Measurement of non-commuting spin components using spin-orbit interaction. 
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We propose a possible experiment aimed at a joint measurement of two non-commuting spin 1/2 
components and analyze its physical meaning. We demonstrate that switching of a strong spin- 
orbit interaction, e.g., in a solid state or a cold-atom system, for a short time interval simulates 
a simultaneous von Neumann measurement of the operators ax and ay. With the spin dynamics 
mapped onto the quantum coordinate-space motion, such an experiment determines averages of 
ax and ay over the duration of the measurement, however short the latter may be. These time 
averages, unlike the instantaneous values of ax 
arbitrary accuracy. 



and cTy, may be evaluated simultaneously to an 
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Recent developments in quantum information and 
technology have brought the quantum measurement the- 
ory (QMT), originally formulated together with the prin- 
ciples of quantum mechanics (see, e.g. [l|) into the re- 
search focus. Continuing progress in experimental tech- 
niques has made it possible to test the QMT as well as 
make new, sometimes surprising, predictions One 
fundamental problem in the QMT is that of joint mea- 
surement of non-commuting variables which, according 
to the uncertainty principle, cannot have well-defined val- 
ues simultaneously. An operational approach to the joint 
measurement of particle's position and momentum was 
roposed in the pioneering work of Arthurs and Kelly ^ 
. Recent attempts to extend it to non-commuting spin 
components can be found in Refs. 0-0. Still, important 
questions concerning the exact nature of the measured 
quantity, the accuracy to which its value can be deter- 
mined, and the back-action a exerted on the measured 
system remain unanswered to this day. The purpose of 
this Letter is to answer these questions, crucial for under- 
standing the nature of any quantum measurement. We 
also suggest an optimal experimental technique for sim- 
ulating a joint von Neumann measurement on a generic 
spin- 1/2 system, of interest in quantum information. For 
the latter we propose the use of modern techniques devel- 
oped for controlling spin-orbit (SO) interactions in solids 
[lo[ and for cold atoms in optical lattices pj]-[l^ . The 
key feature of such systems, currently attracting inter- 
est for both fundamental and applied reasons (for a re- 
view see 0), is entanglement between the translational 
and the spin (pseudospin) degrees of freedom. Gener- 
ated once the SO coupling is switched on, the entangle- 
ment allows the particle play the role of a von Neumann 
pointer. Modulation of the SO coupling strength, includ- 
ing switching it on and off on demand, can be achieved 
for electrons in semiconductor structures by applying ex- 
ternal bias to the metallic gates attached to the system 
S [iol • For cold atoms similar effect can be realized with 
specially designed optical fields pj]-[l^. 
With the above in mind, for a particle of mass M, we 



will consider one of the following Hamiltonians {h = 1): 
H = g{t){p,a^ ±Py(T6) +pV2M, (1) 



where Px and Py are the components of two-dimensional 
momentum, = + Py and the indices of the Pauli 
matrices a, are either j = x^6 = yoT'y = y^6 = x. We 
assume the SO interaction to be switched on for a finite 
period of time, 

g{t) = a for < t < T, g{t) = otherwise, (2) 

where the coupling parameter a varies from 10 cm/s 
for cold atoms to 10^ cm/s for electrons in semiconduc- 
tors. Thus, for < t < T, the operator of particle's 
velocity depends on the orientation of the particle's spin, 
'^x = Px/M -\- acr^, Vy = Py/M ^ aas. Without loss of 
generality, in Eq.([l]) we choose the SO coupling in the 
form g{t){p ' a) [15|. Neglecting the kinetic energy (ex- 
act condition will be given further in the text) we have 
the Schroedinger equation 

idt^{x, y, t) = -ia{dx(Jx + dy(jy)'^{x, y, t), (3) 

with an initial condition: 



(4) 



where ^{x^y^t) and ?7[^"] are two-component spinors. Us- 
ing translational invariance of the Hamiltonian ([1]) we 
rewrite (|3]) as 



^{x,y,t) 



G{x -x\y- y')r]{x', y\ t)dx'dy', (5) 



r){x,y,t) = {x\{y\eM-iHt)\0)\0)r]^' 



(6) 



where, in addition, J rj{x,y^t)dxdy = ?7[^"]. We note 
that Eqs.(|3])-(|6|) are identical to those describing a spin 
coupled to two von Neumann pointers [l| with positions 
X and y^ respectively, which attempt to measure two 
non-commuting projections of the spin simultaneously. 
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FIG. 1: (color online) (a) two possible particle's virtual paths 
in the x^y-plane leading to the same final states (b) the solid- 
line path in the dx-subspace with (ctcc)^ = 1/3, (c) same path 
in the cr^-subspace with (cr^)^ = 0. 



Considering first measurement of a single spin compo- 
nent, say (Jx^ by choosing H = —iadxCJxi offers a use- 
ful insight. Since operators exp(— at^icr^) effect trans- 
lations in the coordinate space, for ?7[^"] not an eigen- 
state of cFx-, ^{x^t) is split into two components which 
travel in opposite directions with speed a. In this way 
one is able to measure cjx to an accuracy determined by 
the coordinate width of the initial pointer state. Since 
(Tx and cFy do not commute, in a simultaneous measure- 
ment the particle which plays the role of the pointer can- 
not acquire a well defined velocity. To study its mo- 
tion, we slice the time interval [0,T] into L subinter- 
vals e = T/L, send L to infinity and apply the Lie- 
Trotter product formula [l6| to write eyi-p[—aT{dx(Jx + 
dyCjy)] = [eyi-p{—a£dx(7x) exp(— aeS^yCr^y)]^. Using the 
spectral representation for each Cartesian component, 
exp(-a£Sicri) = |m).exp(-ma£:Si)i(m| where 

= m\m)- and i = x, one readily sees that a 
pointer undergoes a virtual random walk on a lattice 
x{jx) = jx(^e, y{jy) = jyae, jxjy = ... - 1,0,1,... 
reminiscent of Feynman's checkerboard for a Dirac elec- 
tron . In every time step the particle moves forwards 
or backwards along the x- and ?/-axes. Its final posi- 
tion is determined by the differences, Aux and An^, be- 
tween the numbers of forward and backward steps taken 
in each direction or, more precisely, by the interference 
between all spacial paths sharing the same Aux^ and 
Ariy (see Figdja)). Next we assign values m{l) = ±1, 
/ = 1, . . . , L to (Ja; in each of the subinterval steps, write 
asAux = as Xl^i ^(0 same for ay. We note 

that finding the pointer at a location (x, y) one also de- 
termines time averages of the spin components, {(Jx)t 
and {(Jy)j^^ defined for the spin-space Feynman paths 
(FigslHb) andlHc)), 



/ crx,y{t)dt, 
Jo 



to an accuracy determined by the position spread of the 
initial state G{x^y). We note further that a particle ini- 
tially localized precisely at the origin can advance along 
the X-axis at most by aT provided all L steps are taken 
in the positive x-direction. Since the Hamiltonian in (j3j) 
is invariant under rotations of the coordinate axes, this 
also implies that the particle would never leave the 'al- 
lowed' circle r = (x^+^^)^/^ < i?so, where we introduced 
the SO coupling determined radius R^o = olT. 
To study the distribution of {(Jx)t (<^?/)t detail 
we return to Eqs.(j3]), where we choose G{x^y) to be a 
symmetric Gaussian of a width tq centered at the origin 



G{x,y) 



V2 
Vttto 



exp(-rV^o)- 



(8) 



It is convenient to define four amplitudes U^j,{x^y^T)^ 
/ii^u = 1,2 for a particle initially at the origin and with 
spin initially polarized along {u = 1) or against = 2) 
the z-axis, to be found at t = T at a location (x,^), 
polarized along (/i = 1) or against (/i = 2) the z-axis. 
Performing a Fourier transform of Eq.(|3]) with respect 
to X and we find, in the cylindrical coordinates, that 
, T) is a Hermitian matrix whose elements are 



kdk 



/7ii(r,T) = /722(r,T) = V2^roX 
J exp (-/c^ro/4) cos(i?so^) Jo{kr)^^, 



(9) 



/7i2(r,^,T) = /7*i(r,^,T): 
J exp {-k'^rl/4) sin {Rsok) Ji (kr) — , 

where Jn(^) is the Bessel function of the first kind of 
order n and is the angle the vector (x, y) makes with the 
X-axis. We are interested in an accurate measurement, 
where the maximum shift of the particle (pointer), i?so, 
is much greater than the width of the initial Gaussian 



ro/Rso < 1, 



(10) 



(7) 



in which case main contributions to the integrals in ([9j) 
come from the region where kR^o ^ 1, and the condition 
for neglecting kinetic energy reads T <C r^M. 

To evaluate possible switching times T for electrons in 
semiconductors, we use recent estimates [9|,[iq| ro ~ 10~^ 
cm, a ~ 10^ cm/s, and the effective mass M ^ 10~^^ g, 
thus obtaining T ~ 10~^^s, meaning that rp / R^g cannot 
be less than 0.1. For cold atoms with [III, [HI ^ 10~^ 
cm, a ~ 10 cm/s, and atomic mass of M ^ 10~^^ g, 
we obtain a broader range of durations, 10 ~^ < T < 
10~^ s, suitable for an accurate simulation of a joint von 
Neumann measurement of the two spin components. 

Replacing the Bessel functions by their large argument 
asymptotes, [20| Jn{z) ^ (2/7rz)^/^ cos(2; - nii/2 - 7r/4) 



3 



and neglecting oscillatory terms in the integrand yields 

(11) 



U{r,e,T)^F{r,T) 
F{r,T)= 



1 exp(— z6>) 

ex.p{iO) 1 



exp 



For a small ro, the radial function F{r^T) shown in 
FiglSfa) has a maximum and a minimum close to r = R^^^ 
rapidly decreases for r > i^so, and exhibits a somewhat 
slower decay for r < R^^. This behavior is understood 
by noting first that for ro = the integral in the first of 
Eqs.(|9]) can be calculated exactly in the Cartesian coor- 
dinates {k = {kl + klY^'^), 



f /7~> 7\T/7\ kdk 
/ cos[Rsok) Jo[kr)—— 

Jo 

dkxdky sm{kRso) 



dRs^ 



J — ( 



k 



(12) 

eyi-^iikxX + ikyy). 



The Fourier transform of sm{kz)/k is known [21] to be 
(27r)~-^(z^ — r'^)~^^'^Xz{^) where the characteristic func- 
tion Xz{^) = 1 for r < z and otherwise. For a finite 
ro, convolution (cf. Eq.(j5])) of ([12]) with the Gaussian (j8j) 
yields an alternative highly accurate form for F(r, T) 



F{r,T) 



1 d 



f 

Jo 



-, exp[-(r-rOVrg] 



dr\ (13) 



and, therefore for /7ii(r, T). From Eq. (p!3|) it is readily 
seen that for small ro, the integral in Eq. (p!2]) behaves as 
{Rso — t)~^I^ if i?so — r > ro, peaks at r ^ i?so when the 
center of the Gaussian coincides with the integrable sin- 
gularity of {R^q—t'^~^I^. For r > i?so it rapidly decays as 
the overlap of the Gaussian with the interval [0, i?so] de- 
creases. Accordingly, F(r, T) behaves as —{Rso — r)~^/^ 
for r < Rso^ passes through a zero at r ~ R^o^ and rapidly 
decays outside the allowed circle r < Rso. Thus, with the 
help of (pT]) for the wavefunction (|5|) we find (indices 1 
and 2 are used for the spin projections up and down the 
z-axis, respectively). 



*i(r,^,T) 
*2(r,^,T) 



= F(r,T) 



r][ -\-r]2 ex.p[—zO) 

[in] /.^N . [in] 

r][ ^ ex.p[zO) + r]2 



,(14) 



concentrated in a narrow ring of a radius Rso and a width 
~ ro- From Eq. (p!4|) we obtain the probability to find the 
particle at a location (r, 6>), 



p(r,^,T) = |F(r,T)|2 



1 



(15) 



where ug = (cos 6>, sin 6>) , and a^y^^ = {r]^^^^\aj\r]^^^^) con- 
sists of the spin components for the initial state (see, e.g., 
[ii,[22|-[2l for solid-state realizations). With the Hamil- 
tonian in (|3]) invariant under rotations in the x^-plane. 
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FIG. 2: (color online) a) F(r, T) in Eq.([n]) for tq/Rso = 0.02. 
Its approximation in Eq. (|13p coincides with F(r,T) with the 
graphical accuracy. This small ro/Rso ratio can be realized 
in cold atomic gases; b) Uii(r,T) (solid line) and Ui2{r,0,T) 
(dashed line) for ro/Rso = 0.2 (realizable in semiconductor 
spintronics) calculated with exact Eq.JO]); c) the spin projec- 
tion av as a function of the coordinate. 




FIG. 3: (color online) Probability distribution p{r,0, T) for 
a spin initially polarized along the x-axis, 77^^^^ = (1, 1)^/a/2 
and ro/Rso = 0.01. 



angular dependence in Eq. (p!5]) comes from the asymme- 
try of ?7[^"] . In particular, for a spin whose initial direction 
is normal to the plane i.e., for ?7[^"] = (1,0)"^ or (0,1)"^ 
the distribution is isotropic, p(r, 6>,T) = |F(r, T)p. Fig- 
ure [3] shows p(r, ^, T) for a spin initially directed along 
the X-axis. 

While in an accurate (ideal) measurement position of 
the particle (x, y) correlates with the time averages {ax)rp 
and {(Jy)j^^ its relation to the final spin orientation is 
less direct. As seen from the Feynman path analysis, at 
some point (x,0), the amplitudes to have polarizations 
along and against the x-axis, build up from the paths in 
Figllja) which arrive at x from the left (e.g., dashed line) 
and from the right (e.g., solid line), respectively. Since no 
paths arrive at x = Rso from the right, the spin at that 
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FIG. 4: (color online) Spin configuration for different forms 
of the SO coupling Hamiltonians: (a) pxCTx -^-Pycry, (b) pxCTy + 

Pyax, (c) PxCTx -PyCFy, aud (d) PxCfy -PyCTx- 

point (and elsewhere on the circle shown in FigHJ^a)) is 
always pointing outwards, 

(r, 6>, T) = cos l9, ay (r, 6>, T) = sin l9, (16) 

where a^(r,6>,T) = (^|a^|^)/(|^ip + |^2p). This is 
also true for any initial spin state rj^^^^ which, since the 
matrix in Eq. (pT|) is singular, cannot be reconstructed 
from the spinor in the r.h.s. of Eq. (p!4|) . Although in 
an ideal measurement the mean pointer's velocity, Vi = 
Vi/T, {(Ti)rj. and a^, are simply related, vi = a{(Ji)j. = 
aai, this is no longer true for a less accurate measure- 
ment. For ro ~ i^so, the final spin state ^(x,?/,T) is 
a superposition of all r]{x\y^^T) which fit under the 
Gaussian in Eq.([5]) centered at {x^y) and the correla- 
tion between Vx^y and ax,y may be lost due to interfer- 
ence. Figure [21(c) shows the projection of the final spin 
onto the particle's mean velocity, a^. The narrow range 
of r's where the particle arrives with its spin reversed 
with respect to its mean velocity, ^ — 1, is a result of 
such an interference. The condition for to equal —1 is 
Uii{r^T) = — [/i2(r, 6>, T) exp(z6>) so that the resonance- 
like feature in Fig|2](c) persists for finite values of ro / i^so 
and vanishes for vq/R^q when the two curves in 
FiglS^b) effectively coincide (cf. Eg. ffTT]) ). 



Our analysis is extended to other Hamiltonians in Eq.([T]) 
by replacing the angle in Eq.(pT]) by — ^ or 7r/2 ± ^, as 
appropriate. Figure [4] shows, for ro/Rso 0, the spinor 
field at t = T for each of the four possible cases. 
In summary, we have shown that switching a strong SO 
coupling over a short time T, in a solid state or cold-atom 
system, simulates a simultaneous von Neumann measure- 
ment of two non-commuting spin components. In this 
case, the particle plays the role of a pointer which corre- 
lates its position, (x, y)^ with the time average of the cor- 
responding spin components, {(Jx)t {^y)T evaluated 
along Feynman paths defined for the two spin variables 
(Figllfb) and [Ifc)). There are infinitely many trajec- 
tories which share the same values of {ax)rp and {(Jy)j., 
leading to the same pointer position. An accurate mea- 
surement reveals that the time averages obey the sum 
rule (ax)'^ + {^y)T ~ ^^^^ angular anisotropy of the 
distribution determined by the anisotropy of the initial 
spin state. Importantly, {(Jx)t {^y)T whose values 
can be determined in a generic joint von Neumann mea- 
surement to an arbitrary accuracy, do not represent 'in- 
stantaneous' quantum mechanical expectation values of 
the operators ax and ay and do not reduce to these, no 
matter how short the measurement is. Indeed, the highly 
irregular fractal-like Feynman paths shown in Fig [1] have 
no intrinsic time scale of their own. Thus, even in the 
impulsive limit T ^ 0, aT = const one does not attain 
unique instantaneous values of the two spin components. 
No matter how short T is, all paths (cf. Figdfa)) con- 
tribute to the transition amplitude ([5]). This is one par- 
ticular way of stating that non-commuting quantities ax 
and ay cannot both have well defined values at the same 
time. 
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